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A three-loop radiative neutrino mass model with dark matter 
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We present a model that generates small neutrino masses at three-loop level due to the existence 
of Majorana fermionic dark matter, which is stabilized by a Z 2 symmetry. The model predicts 
that the lightest neutrino is massless. We show a prototypical parameter choice allowed by relevant 
experimental data, which favors the case of normal neutrino mass spectrum and the dark matter 
with m ~ 50 — 135 GeV and a sizable Yukawa coupling. It means that new particles can be searched 
for in future e + e~ collisions. 
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I. INTRODUCTION 

The discovery of very small, but non-zero neutrino masses and the existence of dark matter (DM) in the Universe 
may provide important information to guide us in the search for new physics beyond the standard model (SM). 
In recent years the idea to incorporate both phenomena in a unified framework has received much attention. And 
among the simplest realizations is the inert doublet model m, which generates one-loop neutrino masses with 
the DM being either an extra scalar-doublet or a Majorana fermion whose stability is protected by an exact Z 2 
symmetry. 

Due to the smallness of the neutrino mass scale, a number of models were proposed to generate neutrino masses 
via higher loop processes, especially via 3-loop ones with the loop suppression (g 2 / 167r 2 ) 3 ~ 10 -13 ( g being a 
electroweak-sized coupling) to naturally explain the large hierarchy m„/f ~ 1CV 13 (v being electroweak scale). 
An earlier model [3] advocated by Krauss, Nasri and Trodden (KNT) extends the SM to include two charged 
scalar singlets and a right-handed neutrino. Meanwhile, the model has an additional discrete symmetry, which 
makes neutrino masses be first obtained at the 3-loop level via the new particles with the masses of order of TeV. 
Therefore, this model is phenomenologically interesting, and is well studied in the subsequent literatures EM]. 
Moreover, the generation of 3-loop neutrino masses also appear in the cocktail model m, which adds to the SM 
two scalar singlets (singly and doubly charged) and a scalar doublet. 

In this paper, we present a new model by substituting a scalar triplet with hypercharge Y = 0 for a charged 
scalar singlet in the KNT model. Similarly, due to the additional Z 2 symmetry and the field content of the 
model, Majorana neutrino masses are also first generated at the 3-loop level, and the lightest Z 2 -odd right- 
handed Majorana fermion could be a DM candidate. 

The paper is organized as follows: in Section 2 we describe the model, obtain the neutrino mass matrix, and 
calculate the DM annihilation processes. Various constraints on the model are analyzed numerically in Section 
3. Then conclusions appear in Section 4. 


II. A MODEL FOR NEUTRINO MASSES AND DARK MATTER 

A. The model 


In addition to SM fields, our model includes several right-handed Majorana fermions Nm, a charged SU(2)l 
singlet scalar S~ and a triplet scalar A with hypercharge Y = 0 


A = 


73 A 0 A+ \ 

-73 A °J- 


(1) 


The number of N, f> will be explained below. Moreover, we introduce a Z 2 symmetry under which the new fields 
are all odd, whereas the SM fields are even. Given the symmetry and particle content of the model, the extra 
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FIG. 1: Three-loop diagrams for radiative neutrino masses. 


lagrangian will be 


AC 


1 

2 Tl ' 


(L^A) 2 ] + (D^S-y D»S~ +iN iR jdN lR 
- (^rn Ni Nf R CN iR + g ia Nf R Cl aR S + + h.c.) - V(A,S~,$), 


( 2 ) 


where C is the matrix of the charge conjugation and the covariant derivatives take the forms 

D^A = ^A-i|[W;r a ,A] , (3) 

+ ig'B^S~ . (4) 


Here r a (a = 1,2, 3) is the Pauli matrix. The scalar potential of the new fields and the SM-like doublet <f> looks 
like 


V(A = -g 2 H &<S> + g, 2 s S + S~+g 2 A Ti-[A 2 } + \ 1 (&<S>) 2 + \ 2 {S + S-) 2 

+A 3 (Tr[A 2 ]) 2 + A 4&$)(S + S~) + A 5 Tr[A 2 ](S+S-) 

TAg^^TrlA 2 ] + (A 7 $ t A$S+ + h.c.), (5) 


where 4 = ir 2 &. 

As Z 2 is exact, A has no vacuum expectation value. After electroweak symmetry breaking, for A 7 ^ 0 the 
charged Z 2 -odd scalars A - and S~ will mix 


m 2 (A ,S ) 


2g\ + X 6 v 2 u 2 A 
±v 2 + >fv 2 ) ’ 


( 6 ) 


where v ~ 246 GeV is the vacuum expectation value of <£. They will give rise to two charged mass eigenstates 


/ H 1 A _ ( cos/3 sin/3 A / A A 
y H 2 ) ~ A — sin /3 cos (3 J \ S~ J ' 


( 7 ) 


Now the extra scalars are H 1 , H 2 and A 0 with masses 


mu 1 < m A o 



+ sin 2 /3 rrijj 


< m H2 ■ 


( 8 ) 


B. Neutrino masses 

Explicitly, the lagrangian in Eq. <[2j) breaks lepton number, and can generate a Majorana mass for the left- 
handed neutrinos. However, The Z 2 symmetry strictly forbids the generation of neutrino masses at either 1- or 
2-loop order, and, therefore, the leading contributions to neutrino masses appear at 3-loop level shown in Fig. [l] 
If the model has a single N R , the neutrino mass matrix will predict two vanishing mass eigenvalues like the 
case in Ref. j4] and contradict the neutrino oscillation data [6j. In order to solve the problem, one can add small 
perturbations to the original mass matrix, add more scalars or right-handed Majorana fermions, and so on. In 
this paper, we employ two right-handed fermions N tR (i = 1, 2) with ijin 1 < mjv 2 , which means that the Yukawa 
couplings gi a can be complex and bring about three physical CP violation phases. However, in the following 
discussion, we leave aside the problem of CP violation for simplicity, so gi a takes real number. 
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For the case of mn 2 > m^o m//,, m,y.. rriw, it is appropriate to neglect the complicated contributions of 
Figs. [l| 6 ) and |Tje). Following the method in [7], we obtain the neutrino mass matrix elements arising from the 
remaining Fig.ll[a) 


( FI^ ) atfi ^ \ 9 i a 9 i i3 ^ ct ^ p I i 


i— 1,2 


where Ii is the three-loop integral 
g 4 sin 2 (2/3)mAr i 


Ii = 


r dr 


x {12 [F 2 (r, m 2 Hl , m \ 2 ) - F 1 (r, m 2 Hl , m 2 H2 )]' 


6 (167r 2 ) 3 m| K J 0 r + m 2 
+ [G 2 (r, m 2 Hl , m 2 H2 ) - G 1 (r, m 2 Hi , m 2 H2 )] 2 - F 2 (r, m 2 Hi , m 2 H2 ) [5F 2 (r, m 2 Hl , m 2 H2 ) 
—6-Fi {r, m 2 Hi , m 2 H2 ) - G 2 (r, m 2 Hi , m 2 H2 ) + Gi (r, m 2 Hi , m| 2 )] | . 


Here four integral functions have been introduced 

P 1 j_ 


22 f 1 x(l — x)r + xvi 2 h 

F\ (r, m Hi , m H2 ) = dx In- 2 - ( m Hl m ^ 2 ) > 

Jo m w 

„ , 2 _2 ^ f'j , (l-x)(xr + m 2 l/ ) + xm 2 f 

F 2 (r, m Hl , m H2 ) = dx In- 5 - m H 2 ) > 

Jo TO w 


Gi(r,m 2 Hl ,m 2 H2 ) — 
G 2 (r,m 2 Hl ,m 2 H2 ) = 


r + 


dx x In 


?(1 — x)r + xm 2 Hi 


" l W JO 
r — rriyy + m 


m 


- (m Hl -t rriH 2 ), 


w 


L H i 


m 


- 1 (1 — x)(xr + rof^) + xm 2 H 

dx x In-o- L — (mHi ~* m H 2 ) ■ 


w 


>0 


m 


w 


(9) 


( 10 ) 


The elements of the neutrino Majorana mass matrix M„ can be related to the mass eigenvalues 

M u = UD„U t with D v = Diag(mi, m 2 , m 3 ), 

where U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) leptonic mixing matrix [13] parameterized by 


C12C13 S12C13 Si 3 e 

U = I —S12C23 — Ci 2 S23Sl3e I,S C12C23 — Si2S23Sl3e I<5 


S12S23 — Ci2C 2 3Si3e 


iS 


ifj 

-C12S23 — Sl2C23Sl3e C23C13 



1 0 0 

0 e iai / 2 0 
0 0 e ia2 / 2 


( 11 ) 


( 12 ) 


(13) 


with Cij = cos 6ij , Sij = sin 6ij. 


Although the numerical results depend on the concrete choice of various parameters in the model, the above 
neutrino mass matrix has the following special structure 


aie a 2e 0 


®le ^1 fi ®lr 


M v = a ltl a 2li 0 a 2e a 2/i a 2r 


(14) 


Ulr 0 


0 0 0 


with a,i a = gi a m a \/T l . Therefore, the mass of the lightest neutrino is zero for Det(M„) = 0. 


C. Dark matter 

When Ni is the lightest Z 2 - odd state, it is stable and can be a WIMP dark matter candidate. For mjv 2 mjvu 
we can safely neglect the effect of N 2 on N\ density. The N\ number density get depleted through the annihilation 
process N-\ (p\)N\ (p 2 ) —► 1+ {p 3 )lp (pa) via the f-channel and u-channel exchanges of H^ 2 . The amplitude for this 
process is 


M af 3 = -g* la gip ( - bm _2 + 7-—^-) u{pa)Plu{p 2 )v{p 1 )P r v(p 3 ) 


t — m 2 H t — m~ H 


+ 3 iq 3 i /3 


sin 2 fj 


i2 Hi 


+ 


cos 


l H 2 


u(pi)P L u(pi)v(jp 2 )P R v(jp 3 ), 


(15) 
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where t = ( p\ — p^) 2 and u = ( p\ — pi) 2 are the Mandelstam variables corresponding to the t and u channels, 
respectively. After squaring, summing and averaging over the spin states, the total annihilation cross section in 
the non-relativistic limit is given by 


^^rel — 




487T 


m N 1 V ie\ 


sin 4 /3 (mj Ii + m+) cos 4 /3 ( m% 2 ± m % i ) 

{m 2 Hi + m 2 Ni ) 4 + (m% 2 + m% i ) 4 


2 sin 2 /3 cos 2 f3 ( m 2 Hl m 2 H2 ± m A Ni ) 
(m 2 Hi ± m 2 Ni ) 2 (m 2 H2 ± m 2 Ni ) 2 


(16) 


where v re \ is the relative velocity between the initial particles. Defining crv le \ = a + bv 2 el , we can approximately 


relate the dark matter relic abundance to the a and b variables by 

1.07 x 10 9 GeV 1 x F 


12 Ay h 


M P 


sfgl a + 3(6 - a/A)/x F 


(17) 


where Mp = 1.22 x 10 19 GeV is the Planck scale, and g* = 86.25 is the number of relativistic degrees of freedom 
at the freeze-out temperature x F given by 


x F = In 


5 /45 g MpniN 1 (a ± 6b/x F ) 

4 V ~8 2t^ 


Here g = 2 is the number of degrees of freedom for the Majorana fermion dark matter. 


(18) 


III. EXPERIMENTAL CONSTRAINTS AND NUMERICAL RESULTS 

Firstly, we summarize some relevant experimental data. A global fit to neutrino oscillation data gives 

si 2 = 0.308 ±0.017, 

s 2 — Q 407+O.O33 /q 4c:c+0.039\ 
s 23 — u -^ * -0.023 l u -^ oo -0.031 ) > 

s 2 — n q2S4 +0 ' 0020 ( 0 0240 +0 '° 019> ) 

s 13 — u -UZo z i_o 0019 V U - UZ ^ U -0.0022/ ? 

A m\ x = 7.54+22 x !0~ 5 eV 2 . 

| Am 2 1 = 2.43 ± 0.06 (2.38 ± 0.06) x 10 " 3 eV 2 , 


(19) 


where the values (values in brackets) correspond to mi < m 2 < m 3 (m 3 < mi < m 2 ), he. normal mass spectrum 
(inverted mass spectrum), and Am 2 = m 2 — (m 2 ±m 2 )/ 2 . As mentioned before, the lightest neutrino in the 
model is massless, thus 

mi ~ 0 (4.89 x 10” 2 ) eV, m 2 ~ 8.68 x 10 ~ 3 (4.97 x 10” 2 ) eV, m 3 ~ 4.97 x 10 " 2 (0) eV. (20) 

The measured value of the relic density from WMAP [16] and Planck m is 

nh 2 = 0.1199 ± 0.0027. (21) 

Moreover, the additional H~ and N, t can mediate 1-loop lepton flavour violating (LFV) processes, such as 
l a —■ y 1/31 (a = p,T, /? = e, fj), and the branching ratios are 


Br (l a -± 7 ) — 


3a 

64t tG 2 f 


^ ) 9ia9i0 

i= 1,2 


sin 13 ( m 


m 


H 1 


m 


H 1 


cos 2 (3 f m 
72 24 ( 7772 


m 


h 2 


m 


h 2 s J 


x Br(Z a ^ Ip^a^f 3 ), 


( 22 ) 


where a = e 2 /(Att) is the electromagnetic fine structure constant, G F is the Fermi constant, and the function 
H{x) is given by 


H{x) = 


1 — 6x ± 3a ; 2 ± 2a ; 3 — 6 a ; 2 In 


6(1 — a;) 4 


(23) 
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FIG. 2: The allowable values of mH 1 and mjv 2 satisfying the neutrino oscillation data, LFV constraints and the observed 
DM relic density. 


The current experimental upper bounds of the LFV processes are mm 


Br(/x —>• ey) 
Br(r —► ey) 
Br(r —> /zy) 


< 5.7 x 1(T 13 , 

< 3.3 x 10" 8 , 

< 4.4 x 10 ” 8 . 


(24) 


In addition, our model can generate effective four-lepton contact interactions at the 1-loop level, which can be 
probed in e + e - collisions. Therefore, precise data from LEP will produce limits on the leptophilic dark matter. 
The detail discussions can be found in Ref. mm- 

Now, we illustrate the allowed parameter space in the case of normal neutrino mass spectrum. The relevant 
parameters in the model can be chosen as four particle masses rriH ,, mjVj, a mixing angle [3 and six coupling 
constants gi a . 

In general, the structure of the neutrino mass matrix in our model is inclined to the hierarchy of gi e > gt^ > gi T , 
and the observed relic abundance implies that \dia9ip\ is of order 0(1 —10) for 50 GeV < my, < 200 GeV. 

Consequently, the Yukawa coupling constants could produce the large LFV branching ratios contradicting the 
current data, especially for g —> ey. However, it is interesting that for suitable parameter values the contributions 
of N\ and IV 2 in Eq. (22) can cancel out due to the opposite sign between g\ e gi M and <? 2 e 52 /*• 

In Fig. [2j we keep m/v, = 100 GeV, g\ e = 0.9, and use experimental data in Eq. (19), (21) and (24) to determine 
the allowable values of m/q and mjv 2 according to mjj 2 and cos 2 /3. To guarantee the expression of the neutrino 
mass matrix in Eq. only considering the contribution from Fig. [l](a), we pick mq 2 > 1 TeV, which means 
that mao in Eq. (|8j) is much larger than TOjq, and rriw- In this figure, one can find that the determination 
of mjq weakly depends on mq 2 . For larger cos 2 /3, such as cos 2 /? > 0.5, Ni is heavier than and can not be a 
DM candidate. Meanwhile, for larger m h 2 ■ such as mn 2 > 7 TeV, g\ e gi^ and g 2 e 92 ^ have the same sign, which 
gives rise to unacceptable LFV. 

In fact, six coupling constants gi a are also definite in the case of Fig. [2] Now we explicitly present their 
prototypical values in Fig. [3] Here we assume 


= 5 TeV, m,H 1 = 1.5iriN 1 


cos 2 (3 = 0.25, 


(25) 


and m,v 2 takes a suitable value to realize the cancellation in the decay g —> ey. As for my > 135 GeV, is 
negative, so the cancellation disappears, which gives rise to unacceptable LFV. From the figure, we can give a 
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m Ni [GeV] 

FIG. 3: A prototypical choice of the Yukawa couplings. Note that gi fl and gr; T are multiplied by 1(V 2 and 10 -3 , respectively. 


benchmark in the model. All input parameters are 


m Nl = 100 GeV, m N2 = 350 GeV, m Hl = 150 GeV, m H2 = 5 TeV, 
cos 2 /? = 0.25, g le = 0.909, g lfl = -4.52 x 1CT 3 , g lT = - 1.29 x 1CT 3 , 
g 2e = 1.01, g 2fJ ,= 1.35 xlCT 2 , g 2r = 5.99 x 1CT 4 , 


(26) 


which leads to the neutrino oscillation data in Eq. (19), the DM relic density in Eq. (21) and the following LFV 
results 


Br(/x —► ej) 
Br(r -» ej) 
Br(r —> /ry) 


5.3 x 1(T 14 , 
1.8 x 1(T 12 , 

1.3 x 10~ 16 . 


(27) 


They are consistent with the experimental bounds in Eq. (24). 

Moreover, discussion in the case of inverted neutrino mass spectrum is similar, but the bigger coupling constants 
lead to a much smaller viable parameter space. 


IV. CONCLUSION 

In this paper, we discussed an extension of the SM which includes two right-handed Majorana fermions JV*, 
a charged SU(2)l singlet scalar S~ and a triplet scalar A with hypercharge Y = 0. Due to the additional 
Z 2 symmetry, the Z 2 -odd fermion N-y could be a DM candidate and generate Majorana neutrino masses at the 
3-loop level. Furthermore, the model predicts that the lightest neutrino is massless for the particular structure of 
neutrino mass matrix. We also analyzed the constraints on the model coming from relevant experimental data, 
and presented a prototypical allowed parameter choice, which favors the case of normal neutrino mass spectrum 
and dark matter with to ~ 50 — 135 GeV and a sizable Yukawa coupling constant g\ e . It means that the DM and 
the charged scalar can be searched for in future e + e~ collisions. 

Finally, we did not consider the problem of CP violation in the model. According to recent analyses J22., ' 251 - 
the best fit value of the Dirac CP violation phase is S = 3tt/2. Therefore, the coupling constants gi a can be 
complex, and the model will possess more phenomenology. 
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